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ABSTRACT 

We analyzed the rapid aperiodic X-ray variability of different types of X-ray binaries 
(black hole candidates, atoll sources, the recently discovered millisecond X-ray pulsar, 
and Z sources) at their lowest inferred mass accretion rates. At these accretion rates, the 
power spectra of all sources are dominated by a strong band-limited noise component, 
which follows a power law with an index roughly 1 at high frequencies and breaks 
at a frequency between 0.02 and 32 Hz below which the spectrum is relatively flat. 
Superimposed on this, a broad bump (sometimes a quasi-periodic oscillation) is present 
with a 0.2-67 Hz centroid frequency that varies in good correlation with the frequency 
of the break. The black hole candidates and the low-luminosity neutron star systems 
(including the millisecond X-ray pulsar) have the same relation between the frequency of 
the bump and the frequency of the break. These similar characteristics strongly suggest 
that in all those different types of sources the band-limited noise and the bump are 
produced by the same physical mechanism. This mechanism cannot then depend on 
the presence or absence of either a small magnetosphere or a solid surface, so that it is 
most likely related to an instability in the flow in the accretion disk that modulates the 
accretion rate. 

The Z sources, which are more luminous than the other sources discussed here, follow 
a similar, but slightly shifted correlation between the break frequency and the frequency 
of the bump. The data suggest that the band-limited noise in Z sources is more complex 
than that in the other sources. 



Subject headings: accretion, accretion disks — stars: neutron, black hole — X-rays: 
stars 
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1. Introduction 

At their lowest observed mass accretion rates (M) 
the low-luminosity neutron star low-mass X-ray bina- 
ries (LMXBs; these are usually atoll sources; Hasinger 
& van der Klis 1989) and the black hole candidates 
(BHCs) are remarkably similar with respect to their 
X-ray spectra (e.g. Barret & Vedrenne 1994; see, how- 
ever, Hcindl & Smith 1998) and their rapid aperiodic 
X-ray variability (e.g. van der Klis et al. 1994a, b; 
Bcrger & van der Klis 1998). So far, the only clear 
observable difference is the detection of quasi-periodic 
oscillations (QPOs) between 200 and 1200 Hz in the 
atoll sources (Strohmayer et al. 1996; see van der 
Klis 1998 for a review) and not in the BHCs (e.g. 
Remillard et al. 1998). Apart from these high fre- 
quency phenomena, the power spectra are (almost) 
indistinguishable. At their lowest observed M the 
power spectra of both the atoll sources and the BHCs 
are dominated by a strong band-limited noise compo- 
nent, which follows a power law with index of roughly 
1 at high frequencies and breaks at a certain frequency 
(14,) below which the spectrum is relatively flat. Su- 
perimposed on this noise component, above i^,, a 
broad bump is present. This broad bump is some- 
times called a QPO although it usually does not sat- 
isfy the commonly used criterium (the ratio of the full 
width at half maximum and the centroid frequency 
is less than 0.5) to be called QPO. It was shown 
that the recently discovered (Wijnands & van der Klis 
1998a,b) millisecond X-ray pulsar SAX J1808.4-3658 
was indistinguishable with respect to the X-ray spec- 
trum (Gilfanov et al. 1998; Heindl & Smith 1998) 
and the rapid aperiodic X-ray variability (Wijnands 
& van der Klis 1998c) from the atoll sources. It has 
been suggested that the same physical mechanism is 
responsible for the rapid aperiodic timing behavior in 
these different types of sources (van der Klis 1994a,b), 
and from the properties of the millisecond X-ray pul- 
sar it was concluded that it most likely originated 
outside the magnetosphere (Wijnands & van der Klis 
1998c). 

Due to the large flow of data generated by the 
Rossi X-ray Timing Explorer (RXTE) satellite, de- 
tailed quantitative comparisons between atoll and 
BHC power spectra are only now becoming possible. 
In this Letter, we present a detailed analysis of the 
low-frequency power spectra of atoll sources, BHCs, 
and the millisecond X-ray pulsar. We show that v\, 
is well correlated with the frequency of the bump for 



the three different source types, strongly indicating 
that indeed the same physical mechanism is respon- 
sible in these sources. We also show that the bright- 
est neutron star LMXBs (the Z sources) at their low- 
est observed mass accretion rates have similar power 
spectra, however, they follow a somewhat different 
correlation. 

2. Observations, analysis and results 

We used data obtained with the proportional counter 
array (PCA) onboard RXTE to study the low-frequency 
power spectra of the different types of X-ray binaries. 
When a strong band- limited noise component was de- 
tected with superimposed on it a broad bump or a 
QPO (hereafter usually referred to as QPO), we fitted 
the 0.004-4096 Hz (or 0.004-2048 Hz when the time 
resolution was insufficient to go to 4096 Hz) power 
spectra. Before fitting the power spectra the dead- 
time modified Poisson level (Zhang 1995; Zhang et al. 
1995) was subtracted. The fit function consisted of a 
broken power law (P oc i/~ ab <= l °™. ab °™ ; where abeiow 
and a a bove are the power law index below and above 
i>h, respectively) to represent the band- limited noise, 
and a Lorentzian with centroid frequency vqpo rep- 
resenting the QPO. In the case of Cyg X-l and GX 
339-4 we used a twice broken power law instead of a 
once broken power law to fit the band-limited noise 
adequately. The second break occurred at frequencies 
between 4.6 and 7.3 Hz, well above the frequency of 
the bump in those sources. Above this second break 
the spectrum steepened further with an index of typi- 
cally 1.8. For some atoll sources we had to include an 
extra cut-off power law (P cx i/-0e~ v ' Vcnt - otf ), repre- 
senting a noise component at frequencies above 100 
Hz (see also Wijnands & van der Klis 1998c; Ford & 
van der Klis 1998). Sometimes kHz QPOs were also 
present, which were fitted with Lorentzians. Errors 
on the frequencies were determined using A\ 2 = 1. 

In many of the power spectra, one or both of the 
desired components were not present. Usually these 
power spectra were obtained when the sources were 
in different source states compared to those where 
the strong band-limited noise and the QPO were de- 
tected. In these states QPOs are sometimes present, 
however, the band-limited noise usually is not. It 
turned out that the power spectra which were suitable 
for our purposes were usually, with a few exceptions 
(see Section ^l]), those taken during the lowest ob- 
served inferred mass accretion rates, i.e., during the 
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low states of the BHCs, during times when the atoll 
sources were in their island states, and when the Z 
sources were at the leftmost end of their horizontal 
branches. Typical power spectra of the source types 
are shown in Fig. [l| 

We also searched the literature for publications 
of similar power spectra as observed with RXTE or 
other X-ray (EXOSAT, Ginga) or 7 -ray (CGRO/OSSE) 
satellites. When such power spectra were encoun- 
tered we used the values for v\, and vapo &s they 
are given in those publications. In Table □ all sources 
which were found and their references are listed. Most 
sources we fitted ourselves with the fit function de- 
scribed above. All but two of the sources found in 
the literature were fitted using the same fit function 
we used; in GRO J0422+52 (Grove et al. 1998) and 
IE 1724-3045 (Olive et al. 1998) the band-limited 
noise was fitted with a zero-centered Lorentzian and 
we estimated the break-frequency from the published 
power spectra. Several other sources (e.g. GRS 1737- 
31: Cui et al. 1997c) also showed the band- limited 
noise, however, no evidence for QPOs was found. The 
very complex sources Cir X-l (neutron star system) 
and GRS 1915+105 (BHC) we left out of our analysis. 
Their complex nature, combined with the sometimes 
many QPOs, makes it difficult to correctly identify 
the QPOs and even the break frequency. 

2.1. The sources 

We have found eight BHCs (Table |l|) which exhib- 
ited the band- limited noise and the QPO. We fitted 
RXTE/PCA data from the BHCs Cyg X-l, GX 339- 
4, XTE J1755-324, and GRO J1655-40 (see also Bel- 
loni et al. 1996, Cui et al. 1997a, 1997b, Revnivtsev 
et al. 1998, and Mendez et al. 1998 for more details 
about the data). From the literature we obtained the 
results of the BHCs XTE J 1755-524 (Fox & Lewin 
1998), IE 1740.7-2942 (Smith et al. 1997), GRS 
1758-258 (Smith et al. 1997), and GRO J0422+32 
(Grove et al. 1998). Usually the BHCs were in 
the canonical black hole low state (GX 339-4; XTE 
J1755-324; GRO J0422+32; IE 1740.7-2942; GRS 
1758-258; Cyg X-l) but sometimes they may have 
been in one or two other previously described states 
(intermediate or very high; XTE J1748-288; GRO 
J1655-40; Cyg X-l; see Mendez & van der Klis 1997 
for a recent description of BHC source states). For 
two sources (IE 1740.7-2942 and GRS 1758-258) a 
second harmonic to the QPO was detected and in 
GRO J1655-40 both a genuine QPO and a broad 



bump are present, not harmonically related (see Sec- 
tion |J). 

We have found eight atoll sources (including the 
millisecond X-ray pulsar; Table which exhibited 
both power spectral phenomena. Most of the data 
of these sources we fitted ourselves (see Mendez et 
al. 1997, Wijnands et al. 1998b, and Wijnands & 
van der Klis 1998c for more detailed description of 
the data), but we also used results published in the 
literature (4U 1728-34: Ford & van der Klis 1998; 4U 
1608-52: Yoshida et al. 1993; 4U 1705-44: Berger & 
van der Klis 1998; Ford, van der Klis & Kaaret 1998; 
IE 1724-3045: Olive et al. 1998). All but one of the 
sources were in the island state; 4U 1735-44 was in 
the lower banana branch (Wijnands et al. 1998b). 

Of the Z sources we fitted data of Cyg X-2 (Wij- 
nands et al. 1998a), GX 17+2 (Wijnands et al. 1997), 
GX 5-1 (Wijnands et al. 1998c), GX 340+0 (Jonker 
et al. 1998; see also Fig. [jji), and Sco X-l (van der 
Klis et al. 1996, 1997) when they were at their low- 
est observed inferred mass accretion rate (thus at the 
leftmost end of their horizontal branches). The QPO 
in this case is the fundamental of the well-known hor- 
izontal branch oscillations or HBO. In these sources 
the second harmonic to the HBO could frequently also 
be detected. We excluded GX 349+2 from our analy- 
sis because so far this source has not exhibited HBO. 

2.2. The QPO frequency versus the break fre- 
quency 

In Fig. ||a, i'qpo is plotted versus obtained 
for the BHCs (black), the atoll sources (red), and 
the millisecond X-ray pulsar (blue) . Although intrin- 
sic scatter is present, ^qpo is well correlated with 
Vh for all three types of sources. The principal rea- 
son for the scatter could be the complex structure 
of the QPO. Considerable substructure, usually be- 
low the main peak, is often present (see e.g. Fig. 
|l|a and b). Moreover, in IE 1740.7-2942 and GRS 
1758-258 two QPOs are detected, harmonically re- 
lated to each other (Smith et al. 1997). It is possible 
that in the other sources also higher or sub harmon- 
ics are present, which are incorporated into the single 
Lorentzian used to fit the QPO, resulting in a QPO 
frequency which is slightly shifted with respect to the 
correct value. We did not plot the data of IE 1740.7- 
2942 and GRS 1758-258 in Fig. |a (but they are 
plotted in Fig. ^o) because it is unclear which one 
of the two QPO harmonics in those sources is sim- 
ilar to the one fitted in the other sources. A clue 
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to resolve this might be that the fundamentals are 
in both sources only marginally consistent with the 
above described relation between and i^qpo , while 
the second harmonics are completely consistent with 
it. 

Another source for which the identification of the 
QPO is uncertain is GRO J1655-40 (also only shown 
in Fig. [2^3) . In this source above v\> both a gen- 
uine QPO near 6.5 Hz and a broad bump near 9.5 
Hz are present (see Fig. 2 of Mendez et al. 1998). 
The frequencies of both phenomena are consistent 
with the relation between and i^qpoj although 
the frequency of the genuine QPO sometimes only 
marginally so. Such confusion might also be present 
in some of the other sources which do not clearly show 
both phenomena, although they could be present. 
This again could result in a i^qpo somewhat shifted 
with respect to the correct value. 

Despite the intrinsic scatter, the relation between 

and i'qpo is very remarkable, extending over three 
decades in frequency for both axes, incorporating 
three different source types (BHCs, atoll sources, and 
the millisecond X-ray pulsar). The ratio of i^qpo to 
does not remain constant, but it decreases with 
increasing frequency. 

We investigated whether the Z sources at their low- 
est observed mass accretion rates also would fit this 
relation. The Z sources are plotted in Fig. [2)3 (red). 
Although a correlation between v\> and ^qpo does 
exist here too, it is clear that the Z sources do not 
follow the same relationship as the other sources. 

3. Discussion 

We have presented results on the broad band power 
spectra of different types of X-ray binaries. In BHCs, 
atoll sources, the millisecond X-ray pulsar, and the Z 
sources similar power spectral shapes below approx- 
imately 100 Hz are observed when they are at their 
lowest observed mass accretion rates. Fig. || shows 
that there is a good correlation between the break fre- 
quency of the band-limited noise component and the 
frequency of the QPO that is often observed above 
this break. The same correlation applies to the BHC, 
the atoll sources, and the millisecond X-ray pulsar, 
a slightly different one to the Z sources. Clearly, in 
Fig. ^) the Z sources are shifted to the left or to 
above the data obtained for the other source types. 
Perhaps the QPOs observed in the Z sources are the 
second harmonic instead of the fundamental, which 



would shift the data points of the Z sources exactly 
on top of those of the atoll sources and the BHCs. 
However, this seems to be ruled out by the tight cor- 
relation, covering atoll and Z sources, between these 
QPOs and the kHz QPOs (Psaltis, Belloni, & van der 
Klis 1998), indicating that the QPOs in Z sources are 
the same phenomenon as those observed in the atoll 
sources. This means that the v\, obtained for the Z 
sources is not the one obtained for the atoll sources, 
perhaps due to the fact that the strong band-limited 
noise components are of different origin in the two 
types of sources. 

3.1. BHCs, atoll sources, and the millisecond 
X-ray pulsar 

The good correlation between and i'qpo i n 
BHCs and atoll sources (including the millisecond X- 
ray pulsar) suggest that the band-limited noise and 
the QPO are caused by one and the same physical 
mechanism (see also Wijnands & van der Klis 1998c). 
The similarities between these source types show that 
the presence or absence of a solid surface and a mag- 
netosphere do not affect these rapid X-ray variabil- 
ity components. These components then most likely 
originate somewhere in the accretion disk at a dis- 
tance of at least several tens of kilometers from the 
central compact object outside a possible magneto- 
sphere, which has a radius of approximately 30 km 
for the millisecond X-ray pulsar (Wijnands & van der 
Klis 1998b). However, the large amplitudes of the 
band-limited noise (up to 50% rms) exclude that the 
emission carrying these fluctuations originates this far 
out in the accretion disk, because most of the gravita- 
tional energy of the accretion disk is released closer to 
the compact object. The apparent lack of inclination 
effects excludes, from a statistical point of view, line 
of sight obscuration as a mechanism. A modulation 
of the accretion rate due to instabilities in the flow in 
the region of the disk outside several 10 km from the 
compact object then remains as the most likely mech- 
anism for generating the band-limited noise and the 
QPO. Another conclusion is that the HBOs observed 
in the Z sources, if they are indeed due to the same 
physical mechanism as the QPOs seen in the atoll 
sources (see above), cannot be explained by magne- 
tospheric beat frequency models (Alphar & Shaham 
1985; Lamb et al. 1985). 

The question arise which fundamental properties 
of the sources determine the exact values of i>h and 
fQPO- The effect of the magnetic field of the com- 
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pact object is probably small in the region of the 
disk where the frequencies are determined. The other 
physical parameters which can affect the accretion 
disk are the mass and spin of the central object, and 
the mass accretion rate. For a given source the mass 
and spin do not change and most likely only the mass 
accretion rate determines and i'qpo (both are 
thought to be positively correlated with M; e.g. van 
der Klis 1994b). Usually the relation between the fre- 
quencies and the inferred mass accretion rate is an 
one-to-one relation, however, it was shown that this 
is not the case for the millisecond X-ray pulsar (Wi- 
jnands & van der Klis 1998c). The reason for this is 
as yet unknown, but it could be related to the tran- 
sient nature of this source (see also Wijnands & van 
der Klis 1998c). The difference between the sources 
with similar mass accretion rates would be the mass 
or the spin rate, or both, of the compact object. This 
could explain why on average the BHCs have smaller 
frequencies than the neutron star systems. However, 
considerable overlap between these source types oc- 
cur, indicating that the mass accretion rate differ- 
ences dominate the frequencies. 

3.2. The Z source Sco X-l in more detail 

When examining the power spectra of Sco X-l, an 
extra noise component is present in the frequency re- 
gion between the band-limited noise and the funda- 
mental of the QPO (van der Klis et al. 1997; see also 
Fig. ||). A similar noise component just below the 
QPO can also be observed in GX 17+2 (Homan et al. 
1998). Perhaps this extra noise component is similar 
to the band-limited noise observed in the other types 
of sources. The dominant band-limited noise in the 
Z sources is then something different, and could be 
related to, e.g., the 6-7 Hz QPO seen in Z sources at 
higher mass accretion rates (see also van der Klis et 
al. 1997). 

In order to test the hypothesis that this extra noise 
component in Sco X-l is similar to the band- limited 
noise in the other source types, we fitted the band- 
limited noise in Sco X-l with a Lorentzian with a cen- 
troid frequency near zero Hertz and the extra noise 
component with a broken power law. The resulting 
data points (u^ between 25.5 Hz and 31.9 Hz; pqpo 
between 41.6 and 46.4 Hz) are shown in Fig. as the 
red filled circles. The points are shifted to higher 
and, taking into account the systematic effects intro- 
duced by using a different fit function, arc consistent 
with the same relation as the other sources. At this 



point, the relation of the band- limited noise and the 
extra noise component in the Z sources to the band- 
limited noise in the other sources is not entirely clear. 
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Fig. 1. — Typical power spectra of the millisecond 
X-ray pulsar (SAX J1808. 4-3658; a), an atoll source 
(4U 1728-34; b), a BHC (Cyg X-1; c), and a Z source 
(GX 340+0; d). The deadtime-modificd Poisson level 
(Zhang 1995; Zhang et al. 1995) has been subtracted. 



Zhang, W., Jahoda, K., Swank, J. H., Morgan, E. H., 
& Giles, A. B. 1995, ApJ, 449, 930 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v4.0. 
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Fig. 2. — The frequency of the QPO versus the break 
frequency of the band-limited noise. In a the black 
points are black hole candidates (Cyg X-l: filled cir- 
cles; GX 339-4: filled triangles; XTE J1748-324: 
filled star; XTE J1755-324: filled squares; GRO 
J0422+32: open triangle), the red points are atoll 
sources, (4U 1728-34: plusses; 4U 0614+09: filled 
squares; 4U 1608-52: open triangles; 4U 1735-44: 
filled triangles; 4U 1812-12: cross; 4U 1705-44: open 
diamonds; IE 1724-3045: filled circle), and the blue 
filled circles are the millisecond X-ray pulsar SAX 
J1808.4-3658. In b the black points are the BHCs, 
the atoll sources, and the millisecond pulsar, the red 
points are the Z sources (Cyg X-2: open circles; GX 
17+2: open squares; GX 5—1: open diamonds; GX 
340+0: open triangles; Sco X-l: plusses), and the 
blue points are the sources (GRO J1655-40: filled 
squares; IE 1740.7-2942: filled circles; GRS 1758- 
258: filled triangles) for which two QPOs (sometimes 
harmonically related) are present (both frequencies 
are plotted). The red filled circles at the upper right 
in b are Sco X-l, however, to determine v\> not the 
band-limited noise was used but the extra noise com- 
ponent at ~ 30 Hz (see text). Error bars are smaller 
than the size of the data points. 
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Fig. 3. — Typical Leahy-normalized power spectrum 
of the Z source Sco X-l (see also van der Klis et 
al. 1997), indicating the extra noise component at 
~30 Hz. On the y-axis not the Leahy power but 
the frequency times the Leahy power (in units of 
[(counts s _1 )(rms/mean) 2 )]) has been used in order 
for the features to be more clearly visible in the power 
spectrum (see also Belloni et al. 1997). Clearly visible 
are both the lower frequency and the higher frequency 
kHz QPO and the fundamental and second harmonic 
of the HBO; the extra noise component is the bump 
below the HBO fundamental. The continuum at high 
frequencies is affected by instrumental effects. 
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Table 1 

Sources with band-limited noise and QPO 



Source 


Typc a 


Break frequency 


QPO frequency 6 


Satellite 


R.cfcrcnccs c 






(Hz) 


(Hz) 






Cyg X-l 


BHC 


0.026-2.9 


0.29-10.0 


RXTE 


1, 2, 3, 4, 5 


GX 339-4 


BHC 


0.040; 0.064 


0.32; 0.37 


RXTE 


1. 2 


XTE J1748-288 


BHC 


5.9 


32 


RXTE 


6 


XTE J1755-324 


BHC 


0.2-0.4 


1.9-3.3 


RXTE 


2, 7 


GRO J1655-40 


BHC 


0.20; 1.28 


0.77 or 1.49; 6.48 or 9.48 d 


RXTE 


2, 8 


GRO J0422 + 32 


BHC 


0.03 e 


0.23 


CGRO/OSSE 


9 


IE 1740.7-2942 


BHC 


0.43 


2.0-f 


RXTE 


10 


GRS 1758-258 


BHC 


0.125 


0.394-f 


RXTE 


10 


4U 1728-34 


A 


1.4-25.7 


8.0-49.2 


RXTE 


11 


4U 0614+09 


A 


0.69-13.9 


6.3-35.1 


RXTE 


1, 12 


4U 1608-52 


A 


0.2-1.0 


1.0-6.7 


RXTE; Ginga 


1, 13 


4U 1735-44 


A 


23.5; 31.9 


61.7; 67.0 


RXTE 


1, 14 


4U 1812-12 


A 


0.095 


0.85 


RXTE 


1 


4U 1705-44 


A 


0.35; 1.67 


2.3; 11.3 


EXOSAT; RXTE 


15, 16 


IE 1724-3045 


A 


0.1 e 


0.8 


RXTE 


17 


SAX .T1H0H. 1 :ji>r„s 


P 


0.26-1.6 


2.4-12.0 


RXTE 


2, 18 


Cyg X-2 


Z 


3.0-9.9 


36.1-55.6-f 


RXTE 


1, 19 


GX17+2 


Z 


1.8-9.2 


30. 3-61.l/ 


RXTE 


1, 20 


GX5-1 


z 


2.1-5.6 


22.3-37. o/ 


RXTE 


1, 21 


GX 340+0 


z 


2.1-4.3 


19.2-34.9/ 


RXTE 


1, 22 


Sco X-l 


z 


5.0-8.0 


41. 6-46. 4/ 


RXTE 


1, 23, 24 



a BHC: black hole candidate; A: atoll source; Z: Z source; P: millisecond X-ray pulsar 

^Whcn two data points are present they are given separately; when more are present only the range of observed frequencies 
is given. The archival, proprietary, and public TOO RXTE/PCA results were obtained by fitting the data (sec Section Lj) ; 
the other results were obtained from the literature. Typical relative errors on the QPO frequencies were 0.1%-1% for thiTZ 
sources and 1%-10% for the other sources; typical relative errors on the break frequencies were 2%-20% for the all sources, 
except for Sco X-l they were typically rul%. 

c l: Archival/proprietary RXTE/PCA data; 2: Public TOO RXTE/PCA data; 3: Bclloni ct al. 1996; 4: Cui ct al. 1997a; 
5: Cui ct al. 1997b; 6: Fox &z Lcwin 1998; 7: Rcvnivtscv, Gilfanov, & Churazov 1998; 8: Mendez et al. 1998; 9: Grove ct al. 
1998; 10: Smith ct al. 1997; 11: Ford &c van dcr Klis 1998; 12: Mendez et al. 1997; 13: Yoshida ct al. 1993; 14: Wijnands ct 
al. 1998b; 15: Berger & van dcr Klis 1998; 16: Ford, van dcr Klis, & Kaarct 1998; 17: Olive ct al. 1998; 18: Wijnands & van 
dcr Klis 1998c; 19: Wijnands ct al. 1998a; 20: Wijnands ct al. 1997; 21: Wijnands ct al. 1998c; 22: Jonkcr ct al. 1998; 23: 
van dcr Klis ct al. 1996; 24: van dcr Klis ct al. 1997 

d A QPO and a broad bump present, see text 

jvj-j estimated from the power spectrum published in the literature 

Presence of a second harmonic 
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